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Background. Research suggests that suicidal behaviour is aggregated in families. However, due to methodological 
limitations, including small sample sizes, the strength and pattern of this aggregation remains uncertain. 

Method. We examined the familial clustering of completed suicide in a Swedish total population sample. We linked 
the Cause of Death and Multi-Generation Registers and compared suicide rates among relatives of all 83 951 suicide 
decedents from 1952-2003 with those among relatives of population controls. 

Results. Patterns of familial aggregation of suicide among relatives to suicide decedents suggested genetic influences 
on suicide risk; the risk among full siblings (odds ratio 3.1, 95% confidence interval 2.8-3.5, 50% genetic similarity) 
was higher than that for maternal half-siblings (1.7, 1.1-2.7, 25% genetic similarity), despite similar environmental 
exposure. Further, monozygotic twins (100% genetic similarity) had a higher risk than dizygotic twins (50% genetic 
similarity) and cousins (12.5% genetic similarity) had higher suicide risk than controls. Shared (familial) 
environmental influences were also indicated; siblings to suicide decedents had a higher risk than offspring (both 
50% genetically identical but siblings having a more shared environment, 3.1, 2.8-3.5 v. 2.0, 1.9-2.2), and maternal 
half-siblings had a higher risk than paternal half-siblings (both 50 % genetically identical but the former with a more 
shared environment). Although comparisons of twins and half-siblings had overlapping confidence intervals, they 
were supported by sensitivity analyses, also including suicide attempts. 

Conclusions. Familial clustering of suicide is primarily influenced by genetic and also shared environmental factors. 
The family history of suicide should be considered when assessing suicide risk in clinical settings or designing and 
administering preventive interventions. 
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Introduction 

Suicide and attempted suicide are leading causes of 
death and morbidity worldwide. According to prog- 
noses by the World Health Organization, approxi- 
mately 1.5 million people will die from suicide in 2020 
(Bertolote & Fleischmann, 2005). Suicidal behaviour 
runs in families (Brent et al. 1996 ; Runeson, 1998 ; Brent 
& Mann, 2005) and the familial transmission appears, 
at least partly, to be independent of psychiatric 
disorder (Agerbo et al. 2002; Qin et al. 2002, 2003; 
Runeson & Asberg, 2003 ; Kim et al. 2005 ; McGirr et al. 
2009;Brezo et al. 2010). 

* Address for correspondence : Dr D. Tidemalm, Department of 
Clinical Neuroscience, Karolinska Institutet, Division of Psychiatry, 
St Goran, SE-112 81 Stockholm, Sweden. 

(Email : dag.tidemalm@ki.se) 



Theoretically, familial aggregation of suicides could 
be explained by shared genes, shared environments or 
both. In fact, results from prior family, adoption, twin, 
molecular genetic, geographic and immigrant studies 
do suggest that familial clustering of suicidal acts 
has genetic as well as environmental causes (Turecki, 
2001; Baldessarini & Hennen, 2004; Brent & Mann, 
2005; Voracek & Loibl, 2007, 2008; Wasserman et al. 
2007). However, due to underpowered studies and 
limitations of data available for study, the strength 
and patterns of genetic and environmental impact 
on familial aggregation of suicidal behaviour remain 
uncertain (Baldessarini & Hennen, 2004; Baud, 2005; 
Bondy et al. 2006; Brent & Melhem, 2008; Brezo et al. 
2008, 2010; Currier & Mann, 2008). 

Three large twin studies yielded estimates of the 
genetic contribution to risk of a suicide attempt 
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ranging from 17% (adjusted for heritability of other 
risk factors) to 55% (Statham et al. 1998; Glowinski 
et al. 2001 ; Fu et al. 2002). Further, adoption studies 
(Kety et al. 1968; Schulsinger et al. 1979 ; Wender et al. 
1986) have demonstrated increased risk of completed 
suicide in biological relatives of adopted-away 
probands with suicidal behaviour. These results are 
usually interpreted as evidence of a genetic contri- 
bution to suicide risk (Roy et al. 1997; Brent & Mann, 
2005; Ernst et al. 2009). However, these children 
may have been subjected to detrimental environ- 
mental effects in utero (Mittendorfer-Rutz et al. 2004; 
Magnusson et al. 2005; Riordan et al. 2006) or in 
early childhood prior to adoption (Roy, 2002; Mann 
et al. 2005; Currier & Mann, 2008; Hawton & van 
Heeringen, 2009). Further, a recent large Swedish 
adoption study of familial risk factors and suicide 
suggested a specific genetic contribution to suicide 
risk (von Borczyskowski et al. 2010). 

Pooled data demonstrate a threefold risk of suicidal 
behaviour (suicide attempt or completed suicide) 
in first- and second-degree relatives of persons who 
attempted or died from suicide (Baldessarini & 
Hennen, 2004). However, this and other meta-analyses 
(Turecki, 2001 ; Brent & Mann, 2005) were limited by 
heterogeneous study designs as well as varying 
and usually small sample sizes. To date, a few large 
register-based studies of familial clustering of suicide 
have been published (Agerbo et al. 2002; Qin et al. 
2002, 2003; Runeson & Asberg, 2003; Mittendorfer- 
Rutz et al. 2008 ; von Borczyskowski et al. 2010). Two of 
these reported a twofold risk of suicide in first-degree 
relatives of suicide victims after adjustment for mental 
disorder (Qin et al. 2003; Runeson & Asberg, 2003). 
However, none of these prior studies reported risks 
across different classes of first-degree relatives or in 
second- and third-degree relatives. Importantly, such 
risk estimates would increase our understanding of 
genetic and environmental influences on the familial 
aggregation of completed suicide (Brent, 2010). 

We aimed to examine familial risks of completed 
suicide in the total population of Sweden using data 
from 1952-2003. We hypothesized substantial familial 
aggregation of suicide and that suicide risk would 
primarily vary by degree of genetic relatedness (Segal, 
2009; Pedersen & Fiske, 2010). 

Method 

Subjects and measures 

We created a large database by linking three national 
longitudinal population registers in Sweden using 
the personal identity number as primary key. The 
latter is a unique persistent identifier assigned to the 



individual at birth or immigration and used for health 
and welfare statistics (Ludvigsson et al. 2009). The 
Cause of Death Register, held by the National Board 
of Health and Welfare, includes all deaths among 
persons registered as residents of Sweden at the time 
of death. Established in 1952, the register covers 
>99% of all deaths from 1961 and onwards, including 
those occurring abroad (National Board of Health and 
Welfare, 2010). Each record contains the date of death 
and codes for causes of death, including definite and 
uncertain suicide, in accordance with ICD-6, -7, -8, -9 
and -10. The National Inpatient Register (maintained 
by the National Board of Health and Welfare) covers 
all in-patient health care (including the few existing 
private hospitals). The register has complete national 
coverage for all psychiatric admissions from 1973 
and onwards. Each record contains admission and 
discharge dates, diagnoses and codes for attempted 
suicide (method used and if definite or uncertain) ac- 
cording to ICD-8, -9 and -10. Diagnoses and codes for 
suicide attempts are determined and recorded by the 
physician at discharge. The Multi-Generation Register 
(maintained by Statistics Sweden) contains infor- 
mation about child-parent relationships for all chil- 
dren born in Sweden since 1932 and who were still 
alive in 1961. Swedish-born children and children who 
immigrated to Sweden with at least one parent and 
were granted citizenship before age 18 years are linked 
to their parents. On 31 December 2004, the register 
contained 7 969 645 index persons together with their 
biological and adoptive parents; in total, 11384649 
individuals. From these child-parent dyads, it is 
possible to construct larger pedigrees, including rela- 
tives at increasing distances (genetic and environ- 
mental) from each index person. 

From the Cause of Death Register, we identified 
all Swedish residents who died from a definite suicide 
(ICD-6 [1951-1957]: E970-9; ICD-7 [1958-1968]: 
E971.9-979.9; ICD-8 [1969-1986], ICD-9 [1987-1996]: 
E950-9; ICD10 [1997-]: X60-84) or an uncertain 
suicide (ICD-8, ICD-9: E980-9; ICD10: Y10-34) between 
1 January 1952 and 31 December 2003 (« = 83951). 
Uncertain suicides (not coded in ICD-6 and ICD-7) 
were included to avoid underestimation of suicide 
rates (Neeleman & Wessely, 1997). Further, we ident- 
ified all cases of a definite or uncertain suicide attempt 
from the National Inpatient Register (ICD-8, -9 and 
-10, codes correspond to those for completed suicide) 
between 1 January 1973 and 31 December 2003. 

Statistical analyses 

We used a matched case-control design to estimate 
familial suicide risks for different classes of relatives 
and other close relationships. We studied relatives at 
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Table 1. Relative risks of suicide in relatives of all probands who committed suicide in Sweden during 1952-2003 (n = 83 951) compared 
with relatives of matched controls, in a total population cohort of 11 384 649 individuals 



Relation to proband 



Number 
of dyads 



Concordant 
pairs 



Matched 
OR (95% CI) 



First-degree relatives 
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12 




< 

♦ 






1 (7 8_3 SI 

J.l ^.o — J.OJ 
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Second-degree relatives 
Maternal half-sibling 
Paternal half-sibling 
Nephew or niece 
Grandchild 


1121888 
1482456 
11247095 
12698233 


34 
34 
154 
247 










1.7(1.1-2.7) 

1.2 (0.8-1.9) 
1.6 (1.4-1.9) 

1.3 (1.2-1.5) 




♦ 

♦ 








♦- 






Third-degree relatives 
Cousin 


20542261 


134 










1.5 (1.3-1.8) 


Non-biological relatives 
Spouse/unrelated partner 
Adopted child 


7538548 
206692 


882 
11 






♦- 




2.3 (2.2-2.5) 
1.3 (0.7-2.5) 




♦ 







Twins 

Twin (monozygotic) 14018 

Twin (dizygotic, same gender) 19 664 

Twin (dizygotic, opposite gender) 21 814 



0 10 20 30 40 50 



15.4 (5.4^3.8) 
3.0 (1.0-9.7) 
5.2 (1.9-14.4) 



Relative risks were calculated with conditional logistic regression by comparing suicide rates in each proband-relative dyad 
type with population control-relative dyads matched 1 : 5 by relative category, gender, birth year and time at risk. The graph 
depicts odds ratios (OR) and 95% confidence intervals (CI) for all dyads (♦). Each individual in the population could appear 
multiple times in different categories (e.g. sibling, child, cousin, etc.) depending on family pedigree. 



varying distances from probands who had committed 
suicide. For instance, for each suicide proband, we 
specified all possible case-sibling pairs consisting of 
the proband and each of his/her full siblings (who had 
or had not committed suicide). For each case-sibling 
pair, we randomly selected five control-sibling pairs 
matched to case pairs by gender and birth year. 
The controls were selected through risk set sampling; 
that is, individuals were eligible as controls if they 
were alive and living in Sweden at the time of the case 
person's suicide, regardless of whether they later 
committed suicide. This matching procedure was used 
for the probands' first-, second- and third-degree 
relatives, as well as monozygotic twins, spouses/ 
unrelated partners and adopted children. First-degree 
relatives comprised parents, full siblings (including 
dizygotic twins) and children. Second-degree rela- 
tives were grandparents, uncles /aunts, half-siblings, 
nephews /nieces and grandchildren, while third- 
degree relatives consisted of cousins. To increase stat- 
istical power, suicidal behaviour (completed suicide 
or suicide attempt leading to hospital care) was used 
as a broader suicidal behaviour phenotype for classes 
of relatives with few pairs concordant for completed 
suicide. 

We calculated odds ratios (OR) using conditional 
logistic regression in proc phreg in SAS version 9.2 



(SAS Institute, USA). Since multiple probands were 
possible in each family constellation, we adjusted 
for the non-independence of probands by computing 
corrected (less narrow) confidence intervals (CI) with 
a robust sandwich estimator (covsandwich option in 
phreg). A similar case-control method was previously 
applied to calculate familial risks of schizophrenia, 
bipolar disorder (Lichtenstein et al. 2009) and violent 
offending (Frisell et al. 2011). For twin analyses, 
asymmetric CI for proportions were calculated with 
proc freq in SAS version 9.2 (SAS Institute). 

Results 

This total population sample included almost 8.0 
million individuals of known parentage, which, 
together with their parents, constituted 11.4 million 
unique individuals. A total of 83 951 individuals com- 
mitted suicide during the study period. Table 1 shows 
the relative risks of suicide in biological and non- 
biological relatives of probands who died by suicide 
compared to risks in relatives of matched control 
individuals. First, risk patterns suggested genetic 
influences on suicide risk; the risk for full siblings 
(OR 3.1, 95% CI 2.8-3.5, 50% genetic similarity) was 
significantly higher than that for maternal half-siblings 
(OR 1.7, 95% CI 1.1-2.7, 25% genetic similarity), 
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despite similar environmental exposure. Genetic im- 
pact was also indicated by monozygotic twins (100% 
genetically identical) having a substantially higher risk 
for completed suicide than had dizygotic twins (50% 
genetic similarity), albeit with overlapping CI. Further, 
cousins, with 12.5% genetic similarity and low or 
marginal immediate environmental exposure to the 
suicide proband, had a 50% higher suicide risk than 
controls (OR 1.5, 95% CI 1.3-1.8). 

Second, shared environmental influences on suicide 
risk were also suggested. Full siblings had a signifi- 
cantly higher suicide risk than offspring of suicide 
decedents (both with 50% genetic similarity but sib- 
lings having a more shared environment, OR 3.1, 95% 
CI 2.8-3.5 v. OR 2.0, 95% CI 1.9-2.2). Also, maternal 
half -siblings had a higher risk of suicide than paternal 
half-siblings (OR 1.7, 95% CI 1.1-2.7 v. OR 1.2, 
95% CI 0.8-1.9, both with 50% genetic similarity 
but the former with a more shared environment 
than the latter). However, this comparison had over- 
lapping CI. 

Spouses were at increased risk of suicide (OR 2.3, 
95% CI 2.2-2.5); this relative risk was higher than, 
for example, that for offspring. However, increased 
suicide risk among spouses is not necessarily caused 
by the family environment shared by the couple, but 
may also suggest assortative mating. This may occur 
because of environmental factors such as social class, 
which affect the likelihood to select a certain spouse, or 
because of the studied phenotype, although not 
necessarily being a conscious or rational choice of the 
mating partner. We were unable to further distinguish 
between the relative importance of the mechanisms 
that might have caused this association. Assortative 
mating could potentially bias heritability estimates 
for suicidal behaviour, but does not lead to different 
interpretations of the familial risks in this study 
(Frisell et al. in press). 

Finally, as a sensitivity test, we attempted to obtain 
more precise risk estimates for selected comparisons 
based on Table 1. Hence, assuming that suicide 
attempts requiring in-patient care and completed 
suicide might be expressions of the same liability 
to suicidal behaviour, we calculated corresponding 
relative risks for suicidal behaviour (completed 
suicide or a suicide attempt leading to hospital care). 
The resulting OR were similar to those for completed 
suicide ; however, since substantially more pairs were 
concordant for suicidal behaviour, CI did not overlap 
as much (see Appendix). Thus, with increased power 
we now found comparisons of monozygotic and 
dizygotic (same-gender) twins (OR 12.3, 95% CI 6.5- 
23.2 v. 2.3, 1.1-5.0) statistically support genetic influ- 
ences on familial aggregation of suicidal behaviour. 
Similarly, shared environmental influences on suicide 



risk were confirmed by maternal half-siblings having 
significantly higher suicide risk than paternal half- 
siblings (OR 1.5, 95% CI 1.3-1.7 v. OR 1.1, 95% CI 0.9- 
1.2). Finally, genetic and shared environmental effects 
were both statistically supported by adopted-away 
biological children and adopted (non-biological) chil- 
dren having higher risks than controls (OR 1.6, 95% CI 
1.2-2.0 for both comparisons). 

Discussion 

This is the first total population study providing 
robust estimates of familial suicide risk in relatives at 
varying genetic and environmental distance from each 
other ; including half-siblings, grandchildren, cousins 
and spouses. Importantly, since this investigation was 
four times larger than any previous study, we had the 
power to identify shared environmental effects on 
the risk of completed suicide. In line with our two 
hypotheses, we found evidence for familial aggre- 
gation of completed suicide influenced primarily by 
genetic and also by shared environmental factors. For 
example, we found a twofold increased risk of suicide 
in children and a threefold risk in siblings of suicide 
probands, compared with corresponding relatives 
of controls. These risk estimates are somewhat lower 
than the average size reported among first-degree 
relatives (children and siblings) in a smaller Danish 
study of a rather comparable design (Qin et al. 2003). 

A few limitations should be noted. The study was 
based solely on register data, although with excellent 
coverage and validity >40 years. Second, longitudinal 
registers are subject to ' left censoring ', or missing data 
before the date the register started. However, since we 
matched for birth year and time at risk, such losses 
would be similar for both case and control dyads and 
not affect the estimates of familial aggregation (i.e. the 
relative risks). Third, despite the large sample size, 
some results were probably non-significant due to 
insufficient statistical power to detect aggregation 
of uncommon events, such as completed suicides. 
Nevertheless, sensitivity analyses using a broader 
measure of suicidal behaviour including attempted 
suicides leading to hospitalization generally sup- 
ported the main findings. Fourth, the matching pro- 
cedure did not include socio-economic status. Fifth, 
familial risks were not easily adjusted for presence 
of mental disorder in relatives. Only data on hospital- 
ization due to mental disorder were available. 
However, since mental disorders are often subclinical 
or treated in out-patient care and suicidality generally 
increases the probability of being hospitalized for a 
mental disorder, this could have induced bias. 

The familial aggregation of completed suicide was 
influenced by substantial genetic and also shared 
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environmental factors. But what are the possible 
mechanisms? The shared genetic component could 
involve liability to impulsive aggression, a propensity 
to react with aggression or hostility when frustrated or 
provoked (Turecki, 2001 ; Brent & Mann, 2006; McGirr 
& Turecki, 2007; Brent & Melhem, 2008). This should 
also be reflected in the familial aggregation of violence 
directed towards others, such as violent crime (Frisell 
et al. 2011). Personality traits including neuroticism 
or neurodevelopmental vulnerabilities involving im- 
paired working memory or executive functioning re- 
sulting in poor problem solving could also be involved 
(Baud, 2005; Brent & Melhem, 2008; Mann et al. 
2009). Gene variants related, for example, to the 
serotonergic-, noradrenergic- and dopaminergic neuro- 
transmitter systems appear to be associated with the 
risk of suicidal behaviour (Brezo et al. 2008 ; Currier & 
Mann, 2008; Ernst et al. 2009); interaction patterns 
(gene-gene, gene-development, gene-environment), 
however, are largely unknown (Rujescu et al. 2007; 
Brezo et al. 2008, 2010; Currier & Maim, 2008; Ernst 
et al. 2009; Roy et al. 2009; Wasserman et al. 2010; 
Fergusson et al. 2011). The genetic effect probably 
involves a large number of single genes or alleles, 
similar to what has been suggested for schizophrenia 
and bipolar disorder (Purcell et al. 2009). However, it is 
important to note that such genetic contributions do 
not rule out that well-designed, environment-based 
interventions might reduce suicide risk. 

On the other hand, our findings also suggested 
shared environmental effects on familial clustering of 
completed suicide. This could involve social learning 
of inadequate handling or communication of frus- 
tration or other negative affects, including those 
following bereavement of a spouse or relative 
from suicide (Agerbo, 2005; Brent & Mann, 2006). 
Transmission of adverse rearing environments could 
also be involved in the aggregation of suicide among 
parents and their children (Dinwiddie et al. 2000 ; Roy, 
2002 ; Brent & Melhem, 2008). Direct role modelling or 
imitation of self-destructive behaviours, for instance, 
between siblings, is also possible although previous 
studies suggest that this effect is unlikely to be sub- 
stantial (Statham et al. 1998; Brent & Melhem, 2008; 
Burke et al. 2010). 

Regarding future research, large epidemiological 
studies with good statistical power are clearly needed, 
since suicide is a rare event, even among relatives 
of suicide victims. Various strategies to gain more 
specific knowledge on how gene variants, develop- 
ment and environment are associated with suicidal 
behaviour have been proposed. These include: 

(a) refining the phenotype for suicidal behaviour; 

(b) considering different possible pathways to suicidal 
behaviour. Concerning the former, endophenotypes 



(intermediate phenotypes between genes and a more 
overt phenotype or outcome), such as certain person- 
ality traits or altered neurocognitive function, should 
be characterized for suicidal behaviour (Brezo et al. 
2008; Mann et al. 2009). Concerning the latter, several 
possible aetiological pathways for suicidal behaviour 
have been suggested, for instance, direct effects from 
psychiatric disorders, effects of early adversity, gene- 
environment correlation (selection into adversity), 
interactions between culture and genes and causal 
loops from genes to environment and back again 
(Kendler, 2010), These two overlapping approaches, as 
well as other strategies, are likely to require further 
development of technology (e.g. neuroimaging and 
-physiology) for more fine-grained assessment of bio- 
logical mechanisms (Wasserman et al. 2010). 

To conclude, this study provided strong evidence 
on the familial aetiology of completed suicide, pri- 
marily caused by genetic but also by shared environ- 
mental factors. The results confirm the importance of 
considering the family history of suicide when assess- 
ing suicide risk in clinical practice or when designing 
and administering preventive interventions. 
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Appendix. Relative risks of suicidal behaviour (completed suicide or suicide attempt leading to hospital care) in relatives of all probands who 
committed suicide in Sweden during 1952-2003 (n = 83 951) compared with relatives of matched controls, in a total population cohort of 
11384 649 individuals 



Number Concordant Matched odds 

Relation to proband of dyads pairs ratio (95 % CI) 



First-degree relatives 
Child 

Adopted-away child 
Second-degree relatives 

Maternal half-sibling 

Paternal half-sibling 
Non-biological relatives 

Adopted child 
Twins 

Twin (monozygotic) 

Twin (dizygotic, same-gender) 



13 714 253 
66564 

1 121 888 
1482456 

206692 

14 018 
19 664 



4689 
87 

246 
219 

75 

12 
7 



2.1 (2.0-2.1) 
1.6 (1.2-2.0) 

1.5 (1.3-1.7) 
1.1 (0.9-1.2) 

1.6 (1.2-2.0) 

12.3 (6.5-23.2) 
2.3 (1.1-5.0) 



CI, Confidence interval. 

Relative risks in this table were calculated by the same procedure as those in Table 1. 



